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Abstract 



^^ I The dynamics of symmetry-breaking after a quench is numerically simulated 

^:^ ' on a lattice for the (2+l)-dimensional 0(3) model. In addition to the standard 

K^ ■ sigma-model with temperature-dependent $^-potential the energy functional 

Q . includes a four-derivative current-current coupling to stabilize the size of the 

2^ , emerging extended topological textures. The total winding number can be 

QQ I conserved by constraint. As a model for the chiral phase transition during 

the cooling phase after a hadronic collision this allows to investigate the in- 

^ ■ terference of 'baryon-antibaryon' production with the developing disoriented 

■^ . aligned domains. The growth of angular correlations, condensate, average 

orientation is studied in dependence of texture size, quench rate, symmetry 



o 



C~^ I breaking. The classical dissipative dynamics determines the rate of energy 

^"■^ ' emitted from the relaxing source for each component of the 3-vector field 

f^ ■ which provides a possible signature for domains of Disoriented Chiral Con- 

^ ■ densate. We find that the 'pions' are emitted in two distinct pulses; for 

Mh, sufficiently small lattice size the second one carries the DCC signal, but it 

Qh! is strongly suppressed as compared to simultaneous 'sigma'-meson emission. 

We compare the resulting anomalies in the distributions of DCC pions with 



^ 



X 



probabilities derived within the commonly used coherent state formalism. 
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I. INTRODUCTION 

The observation of clear signatures for the chiral phase transition still remains one of the 
outstanding challenges of hadron physics. There is hope that the extreme energy densities 
located inside the spatial region between receding baryonic slabs shortly after ultrarelativistic 
baryon or heavy-ion collision may provide physical conditions for such a phase transition to 
take place. 

From the theoretical point of view we are dealing with a complex dynamical system, 
strongly interacting quantum fields with many different fermionic and bosonic degrees of 
freedom, far from equilibrium, creating and emitting baryons and mesons from a rapidly 
expanding and cooling spatial volume. Correspondingly, there is a vast theoretical literature 
ranging from cascade models in the colored partonic degrees of freedom to hydrodynamic flow 
models for hadronic currents, which might be applicable at different stages and shed light 
on different aspects during the time evolution of such events. If within the hot fireball chiral 
symmetry indeed is restored then the dynamics of its spontaneous breaking and formation 
of a chiral condensate is one important aspect of the cooling process in which the highly 
excited spatial volume returns back to normal physical vacuum. 

A possible strategy which has proven useful for many other cases of phase transitions in 
microscopically very complex systems is to comprise the essential phenomenology into an 
effective action for an ordering parameter field. For the hadron phenomenology near T = a 
chiral field with spontaneously broken 0(4)-symmetry has become a standard and successful 
tool, so it is natural to try to extend and apply this concept to the temperature region 
where the chiral phase transition is expected. In its most simple version the effective action 
then consists of the (3+l)dimensional 0(4) linear a-model, with a suitable temperature- 
dependent $'^-potential for the spontaneous symmetry breaking. Depending on the ratio 
between relaxation and cooling times, the dynamics of the ordering field then follows the 
cooling quench imposed on the temperature dependence of the potential. As in other cases 
of multi-component field ordering, it is characterized by formation and growth of misaligned 
domains separated by domain walls or other topologically nontrivial structures, depending 
on field and space dimensionality |I[]. Numerical simulations in this classical framework have 
been performed and they confirm the transient formation of domains of false vacuum, in 
which the condensate is approaching its vacuum value while the direction of the aligned 
field still deviates from the surrounding vacuum 0] . The possible influence of quantum- and 
thermal effects has been investigated, mostly in mean- field approximation ||^. It has been 
suggested that pions emitted from these domains of disoriented chiral condensate (DCC) 
carry the misalignment in their isospin multiplicities, and in this way may provide a signature 
of the phase transition [Q. 

A peculiar feature of the chiral 0(4)-field in 3 spatial dimensions is that its winding 
density may be identified with baryon density . This has led to a most interesting and re- 
markably successful model for baryon structure and dynamics at T = 0. The corresponding 
topological arguments |^ have been used 0] to estimate nascent baryon-antibaryon multi- 
plicities in random chiral field configurations at high temperature; of course, the first few 
time steps in an ordering evolution lead to an almost instantaneous reduction of the initially 
large and rapidly fluctuating winding densities. 

Strictly, the topological conservation of the total winding number is limited to a field 



manifold where the point $ = is excluded. In numerical simulations of field evolutions on 
a discrete spatial lattice, however, even topological conservation laws have to be enforced 
by constraint, because arguments based on continuity cannot be applied. Therefore there is 
no principal difference for implementing on a lattice the conservation of winding number in 
the nonlinear or linear realization of the 0(4)-symmetry. This allows to study formation of 
condensate and disoriented domains during the phase transition with simultaneously forming 
baryons plus antibaryons such that the total baryon number is fixed to any desired value. 
It is well known that topological textures play an important role in ordering transitions, 
so it is of special interest to study how the stabilization of extended baryonic structures 
and baryon number conservation may affect the growth laws of DCC domains, and possible 
DCC signatures in the emitted radiation. For this purpose it is necessary to include in 
the effective action terms which stabilize the spatial size of the baryons generated in the 
course of the ordering evolution. Scaling arguments show that an additional four-derivative 
term is sufficient for soliton stabilization in 2- and 3-dimensional 0{n) models. A different 
extension, the inclusion of the anomaly term, has been advocated to increase the probability 
of DCC-formation §. 

The linear cr-model (even with explicite inclusion of quarks) sometimes has been used Q 
to investigate the possibility of DCC-signals by considering spatially averaged fields, only. 
Naturally, within such an approximation the effects we are investigating here, are lost. On 
the other hand, there is no immediate need to explicitly include fermionic degrees of freedom 
in our present approach. 

Before embarking on a full-scale simulation of the (3-|-l)-dimensional 0(4) case it is 
helpful to analyse the relevant questions in the (2+l)-dimensional 0(3) model [|10|. It shares 
all essentials with the higher-dimensional case. Although it represents an efficient tool for 



the description of 2-dimensional spin systems in its own right W^, we here discuss it mainly 
in view of its hadron physics generalization. That is, we denote the order parameter field as 
'chiral' meson field, we call the fluctuations 'pions' with internal 'isospin' components, and a- 
mesons, which aquire mass through spontaneous breaking of the chiral 0(3)-symmetry, and 
we call the domains with small values of the condensate and unbroken symmetry 'bags', in 
which the 'baryonic' winding density is confined. Although there is no compelling connection 
to fix the parameters of the model from its higher-dimensional analogue, we choose them 
in correspondence to simulate relative sizes and masses for baryons and mesons. Various 
features of the static solutions of this model, and their formation in relaxation processes after 



a quench have been studied in iT^JTH] . In the present investigation we discuss in more detail 



the different stages which characterize the long-time ordering evolution under dissipative 
dynamics, the dependence of the relevant times on model parameters, different types of 
symmetry breaking, and quench rates. Our main aim, however, is an attempt to extract 
from the dynamics of the evolution the possibility of a DCC-signal. Generally, the emission of 
radiation from a moving source is determined by quantum-field-theoretic amplitudes which 
require a separation of the time-dependent field into a classically moving part and quantum 
fluctuations. The result is commonly put into the statement that the 'intensity of pion 
radiation is proportional to the square of the classical pion field' [Q . This argument is the 
basis of the interest in DCC-pions. We analyse the validity of this statement under the 
assumption of dissipative dynamics, i.e. assuming that the energy loss of the relaxing and 
ordering classical configurations is carried away by propagating fluctuations. Under this 



assumption we obtain the timing and strength and isospin distribution of DCC-pions from 
purely classical considerations. 

This general dynamical scheme is outlined in sect. 2, while the specifics of the model 
and the observables of interest are defined in sect. 3. The essential features of configurations 
during the course of the transition, the growth of angular correlations, decrease of numbers of 
textures, saturation of condensate and average orientation, formation of baryonic clusters, 
are presented in chpts.4,5, and 6, in dependence on parameters of the model, symmetry 
breaking, and cooling rates. Finally, in sect. 7, we analyse the rates of energy loss with respect 
to possible DCC pion pulses, and extract the isospin distributions from large ensembles of 
events. Conclusions are drawn in sect. 8. 



II. NONEQUILIBRIUM FIELD DYNAMICS AND DCC SIGNALS 

We consider an effective action for an 0(3) vector field $ 

5[#] = / (T[#] - W[$]) d^xdt (1) 



where the kinetic energy density T[#] comprises all terms containing time- derivatives of $. 
We consider the evolution of field configurations, which at some initial time deviate from 
the global equilibrium configuration (the 'true vacuum') within some finite spatial region V. 
The equations of motion 6S/6^ = which govern the classical evolution of the initial non- 
equilibrium configuration describe both, the approach towards the vacuum configuration in 
the interior of that spatial region and the propagation of outgoing (distorted) waves into the 
(initially undisturbed) surrounding vacuum. Although, of course, the classical equations of 
motion conserve the total energy E = J {T + U) (Px, the outgoing waves carry away energy 
from the interior of the spatial region in which the total energy E = T + U initially was 
located. 

In quantum field theory the quantization of the propagating fiuctuations describes mul- 
tiparticle production of pions emitted from the relaxing field configurations. In a separation 

^{x,t) = ^,i{x,t) + 6^ix,t) (2) 

where $c/ comprises all of the large amplitude motion of $, and 5$ contains only small 
fiuctuations around ^d, expansion of S[^] to second order in 5$ 



6^S 



lm^)\i^J''^^'^'^''^''' 



+ / / ^^(" ' ^ \#(.-,t0^#(^-,f0 '[^-'i^^^"^' t")d^x'dt'd^x"dt" + 0{6<^^) (3) 

and variation with respect to 5$ provides the classical expression for the fiuctuating part 
5^{x,t) = 5^^\x,t) + jg{x,t,x',t')j^^l^^/x'dt'. (4) 



The homogeneous part 5$'°% after quantization, represents the scattering of field quanta off 
the classical configuration $cz (i-e., these are on-shell distorted waves). The Green's function 
Q{x,t,x',t') of the operator S'^S/{6^6^)\r^ -, relates emitted (i.e on-shell) radiation to the 
source term SS/6^\,^ ,. There is, however, no unique way to separate the propagating 
fluctuations 5$ from a more or less smoothly evolving classical conflguration $c« because 
both, $c« and 5$, are integral parts of one and the same evolving order-parameter fleld 
$. Conclusions drawn from one part only are subject to the arbitrariness of the chosen 
separation. In any case, in an equation of motion that separately describes the evolution 
of $c« we expect a dissipative term to account for the loss of energy through the emitted 
radiation: 

U,, = 5S/5^\^^^^y (5) 

For sufficiently small values of the relaxation constant r the dissipative term dominates the 
time evolution of $c/, propagating parts get damped away and we can replace (H) by the 
Time-Dependent Ginzburg-Landau (TDGL) equation 

l^ci = ~SU/S^\^^^^^. (6) 

The potential energy functional U contains no time derivatives of $, therefore $cz has no 
propagating parts, i.e. it does not pick up field momentum, and it provides a very slowly 
moving adiabatically evolving classical background field. 

Back-reaction of the fiuctuations on the moving classical field can re-enter in different 
ways: as loop corrections they can lead to renormalization of the coupling constants in the 
functional U which thus may become time- or temperature-dependent and drive the spon- 
taneous symmetry breaking; an appropriate noise term to be added on the right-hand side 
of d^) accounts for stochastic interactions of the fluctuations with $ci at flnite tempera- 
ture. Ideally, through that procedure one hopes to achieve a separation (H) such that the 
fluctuating parts 5$ contain no exponentially increasing amplitudes during the relaxation 
process. 

In the following we assume that the coupling constants which appear in the TDGL- 
equations have their renormalized physical values with an appropriate time- or temperature 
dependence and we will restrict the discussion to the consequences of those equations. Evo- 
lutions via d^) of individual initial $c/ conflgurations into flnal quasi-stable conflgurations 
are denoted as 'events', and we shall consider ensembles of events sampled with different 
initial conflgurations. The individual events show a variety of physically interesting features, 
textures, walls, disoriented domains, growth laws, which form the main body of the subse- 
quent discussions. The lack of second-order time derivatives in (j^) means that also localized 
stable structures created in $ci will not propagate but only drift in response to their mutual 
interaction. This is as expected because uniform motion of textures represents zero modes 
which may still be contained in 5$. With (j^) the rate of decrease of the energy Ud = U[^ci] 
stored in the classical configuration ^ci then is 

ilciit) = -^ J \^d{x,t)\'d'x. (7) 



If particle emission from the evolving source is the only mechanism for energy nonconser- 
vation this rate of energy loss has to match the energy flow carried away by the emitted 
radiation. This requirement determines the relaxation constant r in (|^) which fixes the time 
scale for the approach of $c/ towards equilibrium. (Of course, external interference with 
the system like Bjorken rescaling of the spatial metric may provide other sources of energy 
nonconservation. ) 

In quantum field theory, for field fluctuations coupled to a time-dependent classical source 
s{x,t), the probability density for the emission of Uj, field quanta with momentum k and 
mass m is given by 

nfc= |s(fc,Wfc)|V(2^fc) (8) 

where s{k,uJk) is the energy-momentum-space transform of the source, expanded in terms 
of the complete set of distorted scattering states with asymptotic momentum k, and taken 
at the on-shell energy u = Uk = yW+m?. 

In our present case according to (^) the time- dependent source is given by 
SS/S^ I r^ ,j, ^N, . Consistent with our assumptions about the slow motion of the classical 
field $ci we again may approximate here 6S/6^ by —6U/6^. With ^^i obeying the TDGL 
equations (|^) the source field then is just the time-derivative of the classical field ^ci{x, t)/T 
and we obtain for the average number of particles emitted with momentum k and isospin 
component i (omitting the index cl) 

|^.(fc,^.)/rP _ (9) 

2uOk 

With particle emission as the only mechanism for energy nonconservation the total energy 
/S.E = J uJkJ2i'>T'i,kd'^k carried away by radiation during the ordering process must balance 
the loss of energy stored in the classical field. With (|^ this requirement leads to 

^ J \^,i{k,LUk)\^d^k = r I \^,i{x,t)\^d^xdt. (10) 

Note that on the left-hand side of this equation only on-shell frequencies appear, while the 
right-hand side (when written in Fourier space) contains an independent integral over all 
Fourier components of the moving classical field. Of course, together with (^, this relation 
is an implicit definition of the relaxation constant r. 

For the adiabatically evolving process we consider the time t as parameter such that we 
may define 

e,(t) = - U,{x,tfd'x = —i^Kt)) (11) 

T J T 

as the energy carried away per time interval by particles emitted with field orientation i. 
Here " ( )" denotes the lattice average over an L x L lattice. Due to the slow motion of the 
source we expect the particles to be mainly emitted with low energies, i.e. with low momenta 
/c ~ and Ui ~ rrii. In our present context we choose the 0(3)-symmetry to be spontaneously 
broken in 3-direction (selected by the surrounding true vacuum boundary condition or by 
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small explicit symmetry breaking, or by a small bias in the initial configuration). Then 
the i=l,2 ('isospin') components of the 0(3) chiral field constitute two 'pionic' fluctuational 
fields S^i = TTj with small mass mi^2 = ""^tt due to explicit symmetry breaking, while the 
cr-fluct nation 6^3 acquires a large mass m„ due to the spontaneously broken symmetry. So 
we expect the emission rate for pions with isospin component i = 1,2 

nm = ^ (12) 

During the early stages of the ordering evolution the spatial averages (<l>^(t)) will be similar 
for both isospin directions, but at late times with the formation of larger disoriented domains 
they might differ appreciably. The relative abundancies 

m = ^"^f ,,, (13) 

ni{t) +n2{t) 

are free of unknown constants and could serve as DCC signature, if it is possible to separate 
in each individual event the small number of late time pions from the background of those 
produced during earlier stages, from decaying a's, and from other sources. 



III. THE MODEL 

The features which are of physical interest in the chiral phase transition for a cooling 
hadronic gas in the 3+1 dimensional 0(4) model can be explored and demonstrated in the 
lower-dimensional 0(3) model. It is defined in terms of the dimensionless 3-component field 
$ = $$ with $■$ = !, with the following lagrangian density in 2 + 1 dimensions 

C = ^d,^d^^ -^{'^'- '^f'^' + - ^^' P^P' + ^'*3. (14) 

In addition to the usual linear sigma-model with standard $^ potential for the modulus field 
$ this lagrangian contains a four-derivative current- current coupling p^p^^ for the conserved 
topological current 

pM = ^gA'^p^ . (5^$ X a„$), (15) 

Stt 

which satisfies (9^p^ = 0. In the present 0(3) model this nonlinear coupling comprises the 
four-derivative 'Skyrme'-term and the sixth-order current-current coupling of effective chiral 
0(4)-models. 

The independent strengths of the $'^- and the Skyrme couplings are written in terms of 
one common dimensionless parameter A, and a length I which may be absorbed into the 
space-time coordinates. So, for A fixed, d sets the scale for the spatial extent of localized 
static solutions, in lattice implementations relative to the size of the basic lattice unit cell. 
The static solutions have been investigated and described in ||T2|. For continuous coordinates 



(or, on a lattice, for £ 3> 1 ) the total energy (apart from the explicitly symmetry-breaking 
m^ contribution) of such localized solutions is independent of £ (as long as £ is small as 
compared to the total lattice size L). In the following we put the lattice constant to unity. 



If we consider the field $ as the basic order parameter to reflect the chiral phase transition 
as function of the temperature T, then the factor /^ which multiplies the square of the field 
vector in (JT^) drives the spontaneous symmetry breaking. As the temperature goes to zero 
it is supposed to converge towards its vacuum value /^(T — ;> 0) — ^ 1, and it should become 
very small /^(T — > Tc) ^ 1 as T approaches the critical temperature Tc. For very high 
temperature p{T > Tc) is supposed to become negative. 

A small symmetry-breaking term with strength m^ has also been added to the otherwise 
0(3)-symmetric lagrangian ([1^ . This explicit symmetry breaking is applied in 3-direction 
of the internal space. If unequal zero it singles out the 3-component $3 of the $-field as 
that field which aquires the nonvanishing condensate at T = 0, and we will denote this 
component as "cr"-field, while the 1- and 2-components represent the "pionic" (Goldstone-) 
degrees of freedom of the remaining unbroken symmetry. 

It is convenient to have also for nonvanishing rn^ the minimum of the (T = 0)-potential 
located at the vacuum field configuration $ = (0,0, 1), such that the length $ of the field 
vectors approaches unity for T = also for nonzero symmetry breaking. This is achieved 
by subtracting m^£^/A from /^ in the coefficient of the second order term in the potential 
\/($) in (|1D 

f = f{T) - ^. (16) 

This form shows directly how the explicit symmetry breaking supplies the pionic Goldstone 
bosons with mass m^. For T < Tc the temperature dependence of m"^ can be ignored. On 
the other hand, the mass of the a-field is given by 

ml = 2^j\T)+ml (17) 

which is mainly determined by the temperature dependence of P{T). 

Evolution of $c« proceeds through the TDGL equations of motion as they result from 
the potential part of ([T^), (i.e. also from the current-current coupling only the density term 
Pq is kept). As no other time- derivative terms appear on the right-hand side, the constant r 
which multiplies the left-hand side of (^) can be set to one, i.e. in the following the relaxation 
time T serves as time unit. Omitting for convenience the index 'cV we have 

^$. = A$, -j,{^'- f) $. - ^^' ^ + <^^^ + e. (18) 

As far as the dissipative term originates from elimination of fluctuational modes which 
simulate a heat bath for the evolving field configuration, addition of fiuctuational noise ^ is 
required with a strength dictated by the dissipation-fiuctuation theorem. Such a term has 
also been added in (|18[). 

In lattice implementation, we impose periodic boundary conditions on the field configura- 
tions. This implies compactification of coordinate space to a torus S^ ^ S"^ . These boundary 
conditions preserve the global 0(3) symmetry of the model (for vanishing m^). Having 
singled out the 3-direction as the cr-direction in which the T = condensate is formed we 
can impose the stronger vacuum boundary conditions $ = (0, 0, 1) for all points on the lat- 
tice boundary throughout the whole evolution. This implies compactification of coordinate 
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space to the two-sphere 5*^. In both cases the winding density po satisfies J pod'^x = B with 
integer winding number B. Of course, these vacuum boundary conditions represent another 
way of exphcit symmetry breaking which eventually drives the evolution of the condensate 
into the 3-direction. It is, however, a very weak form which is only active at the surface of 
the spatial region, unlike a nonvanishing m^-term which is locally active at every point in 
space. 

With the equations of motion (|T^) being purely first order in time derivatives it is suffi- 
cient to specify for initial conditions at t = only the field configurations themselves. As the 
time-dependent fiuctuations are not part of $ci the high-temperature initial configuration 
is, (apart from a small optional bias in 3-direction due to an explicit symmetry breaker), 
$c/ = in the interior of the spatial region in which the hot chiral field is located and 
takes the true vacuum values on its boundary. The first step in the time evolution is there- 
fore governed by the stochastic force $, alone. But this first step will immediately create a 
configuration $c/ which is of stochastic nature itself. Therefore, choosing each of the three 
cartesian field components at each interior point of the lattice from a random Gaussian 
deviate, and fixing <^ci = (0, 0, 1) on the boundary of the lattice, provides convenient initial 
configurations with initial spatial correlation length less than one lattice unit. From these 
random configurations, those with the desired total winding number B are selected and 
propagated in time. During the evolution of each individual event the fiuctuational noise $ 
at every time step and every lattice point is picked from an appropriate Gaussian deviate. A 
large number of events generated in this way constitute the statistical ensemble from which 
ensemble averages can be obtained at some point in time during the evolution. 

For configurations where the orientation of the vector field $ shows strong local varia- 
tions (as it happens near localized textures or throughout random initial distributions) the 
implementation of the winding density po on the lattice in the differential form (|I3|) is not 
sufficiently accurate to add up to integer total winding number B. Therefore po and its 
functional derivatives have to be evaluated according to the geometrical meaning of po as 
surface area of spherical triangles cut out on the two-sphere S'^ by the image of elementary 



triangles in coordinate space. Details of this procedure have been discussed in |]13| . It allows 
to assign a value p{i,j) of the winding density po to each point {i,j) of the square lattice (or, 
more precisely, to each lattice cell with lowerleft corner {i,j)) such that the total winding 
number 

B=J:pi^,J) (19) 

i,j=0 

summed over the whole lattice is integer and initial configurations can be selected with some 
desired integer value of B. We also define 

D= j:\piz, j)\ (20) 

i,j=0 

by summing up the absolute values of the local winding densities. Of course, for random or 
slowly varying smooth configurations D generally is not an integer. However, if a configu- 
ration describes a distribution of localized textures (and antitextures) which are sufficiently 
well separated from each other, then D is close to an integer and counts the number of these 
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textures (plus antitextures). In that case we can define the numbers N+,N- of "particles" 
and "antiparticles" through 

B = N+-N_ D = N+ + N^ . (21) 

In the following we will briefly call D the "particle number", even if it is not integer and 
comprises textures and antitextures. 

In a lattice implementation the local update of $(z, j) at some timestep at some lattice 
point {i,j) will occasionally lead to a discrete jump in the total winding number B. For 
well-developed localized structures this corresponds to unwinding textures or antitextures 
independently, such that B decreases or increases by one or more units. This will eventu- 
ally happen even for implementations of the nonlinear 0(3)-model where the length $ is 
constrained to unity and B is topologically conserved, because the topological arguments 
based on continuity do not apply to the discrete lattice configurations. It can be prevented 
by an (optional) B filter which rejects such a S-violating local field update at that lattice 
point and timestep. This eliminates all independent unwinding processes. Only simulta- 
neous annihilation of texture and antitexture in the same time step remains possible, and, 
as the update proceeds locally at each lattice vertex it can happen only if texture and an- 
titexture overlap. This 5-conserving evolution is characteristic for the nontrivial topology 
of the nonlinear 0(3)-model and in this way can be implemented as an optional constraint 
also into the linear 0(3)-model. In the continuum limit it implies that a vanishing modulus 
$ = of the field vector is excluded. 



IV. ROLL-DOWN AND DOMAIN GROWTH AFTER SUDDEN QUENCH 

In order to get a feeling for a reasonable choice of the two essential parameters i and A in 
(|1^) we study at first the simple case of a sudden quench without any symmetry breaking. 



neither explicit nor via boundary condition, the so-called critical quench. That is, we put 
m^ = and impose only periodicity for otherwise arbitrary field vectors at the boundary of 
a sufficiently large L x L lattice. The sudden quench is realized by exposing the initially (at 
t = 0) random configuration to the T = 0- value of /^(T = 0) = 1 for all times t > 0. For 
this zero-temperature evolution the noise term ^ is also not present. As an example fig.|l| 
shows one individual evolution of a configuration with fixed winding number i? = 0, on a 
lattice with L = 150, with i = 1 and A = 2. 

As significant observables we show the time-dependence of total energy E, particle num- 
ber D, field modulus ($) averaged over the whole lattice, and the angular correlation length 
Rd- The latter is defined as the (interpolated) half- maximum distance R1/2 of the equal-time 
angular correlation function 

c{R)= E E M^,j)-Mk,i)/ E E 1 (22) 

i,j=Ok,l=(} i,j=Ok,l=(} 



(with the A;, / sum restricted such that the distance r = J{k — iy + (/ — jY between lattice 
vertices (i, j) and (A;,/) lies inside bins of unit size around fixed positive integers R). Note 
that the definition (|22| ) contains only the unit vectors $, and not the full field vectors $, i.e. 
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FIG. 1. Total energy E, particle number D, angular correlation length Rd, and the average 
length of the field vector ($) as functions of time after a sudden quench (for i = 2,\ = 1). The 
total winding number is i? = by constraint throughout the evolution. The dashed line indicates 
the classical growth law Ro ~ t^-^. 



C{R) measures only the angular correlations in a given configuration. As discussed in |T3| 
this definition of Rd serves well to provide an average radius of oriented domains. 

Fig.0 follows the evolution over a very long time (up to t ~ 10^ relaxation-time units) 
and clearly shows the different phases of the ordering process. Up to t ~ 10 the average 
($) stays (or becomes) small (~ 0.1), while the particle number D decreases by almost one 
order of magnitude due to local angular reorientation, also reflected in Rd which increases 
approximately as t°'^. Once the oriented domains have reached sizes which cover several 
lattice units the roll-down phase sets in where $ approaches unity within these domains 
(while remaining small at their boundaries where strong angular variations prevail). This 
roll-down phase takes a few tens of time units and it is accompanied by a significant energy 
reduction and a more rapid increase of Rd, such that by the end of this roll-down the 
oriented domains cover areas with typical radii of 10 lattice units. During the end of, and 
after this phase the isolated textures assume their stable shapes with the appropriate bag and 
angular field profiles. This is accompanied by further angular ordering of the topologically 
trivial field configurations in which these textures are embedded, which proceeds by moving 
boundaries between differently oriented domains. It leads to further decrease of E and 
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increase of Rd, which follows very nicely the classical t°'^ Allen-Cahn law [Q. This phase 



lasts up to several hundreds of time units, without further reduction of the particle number 
D. Due to the existence of the remaining stable bags the average $ does not reach the 
vacuum limit ($) = 1 and due to the floating angular textures the increase in Rd begins to 
saturate after t ~ 10^. On a timescale which is still larger by another order of magnitude the 
stable isolated textures may meet and annihilate as can be observed in fig.0 near t ~ 10^. 

Fig.0 is obtained for B constrained to B = and for scale i = 2, such that resulting bag 
structures extend only over a few lattice cells. The long-time behaviour after the roll down 
phase is higly sensitive to the size of the developing stable textures, i.e. to the fraction of 
lattice space occupied by the bags because, evidently, overlapping textures quickly combine 
or annihilate. Their size is governed by i, and in the range up to £ < 5 the structures created 
are reasonably isolated. We shall see later that for the point of view of DCC signatures the 
long-time behaviour after the roll-down is irrelevant. Defining the typical roll-down time 
tfiD as that point in time when the lattice averaged ($) has reached a value of 0.8, we plot in 
fig.0 tRD as function of i and A, respectively, for 15 events each. For B constrained to B = 
the average values follow power laws Ird ~ ^^'^ and tno ~ A^^'^ with good accuracy. These 
numerical results show that the roll-down time t^D scales as (£^/A)^'^ which means that it 
is essentially determined by the potential term ($^ — 1) alone. It is almost unaffected by 
the fourth-order derivative term which scales with the product (^^A), while the second-order 
gradient terms in U cause only a mild increase of the scaling exponent to the observed value 
of 1.2. 




2 3 4 5 

FIG. 2. The roll-down time i/ju for different values of the scale i (with coupling constant 
A = 2), and different values of A (with scale i = 1), respectively. For each set of parameters 
15 events {B = 0) on a 120 x 120 grid are sampled. The straight lines indicate the power laws 
tRD = 2.6 £^'^ and Ird = 6A"^'^, respectively. 



We interpret the stable textures as baryons floating in the effective meson pion field. 



12 



Fixing their spatial extent in terms of lattice units through the choice of i provides the 
physical scale for the lattice constant. We shall typically choose £ = 2 such that baryons 
carrying one unit of winding number (=baryon number) acquire a radius for their topological 
density of about two lattice constants. With a typical nucleon radius of 0.6 - 1.0 fm, for -^ = 2, 
the lattice constant then corresponds to ~ 0.3 — 0.5 fm. For an accurate description of the 
texture and bag profiles of individual baryons such a small value £ = 2 is not really sufficient, 
as has been discussed in detail in the investigation of static solutions in []T3[. But within the 
present context this aspect is not essential. For physical vacuum boundary conditions this 
choice also suggests the typical lattice size L, which should be sufficiently large to cover the 
spatial extent of the high-temperature area after a typical hadronic (heavy-ion) collision. So 
values of 20 < L < 100 seem appropriate, and they are technically easy to handle. 

Through (0) the value of A is directly tied to the cr-mass. For 1 = 2 the choice A = 2 
leads to ml. = 1 (in inverse lattice units), which appears quite reasonable, being of the order 
of a typical baryon mass. In the discussion in [|TB[ we have seen that for A = 2 the stable 
textures are confined into the interior of well-developed bag structures which we consider an 
essential feature of physical baryons. For all these reasons we shall in the following present 
results mostly for £ = 2 and A = 2. 

In order to obtain the roll-down times for that case we extract from fig.^ the values of 
tjiD for this two-dimensional 0(3)-model 

tRD^G fyj . (23) 

For i = 2, X = 2 this is about 14 relaxation time units. 

To guide the eye for the increase of the correlation length Rd we have included in fig.|I| 
the power law R = 1.3 t°'^ as straight dashed line. Although there are typical deviations 
prior and during the roll-down phase it is remarkable that this classical scaling expectation 
serves well to describe the growth of oriented domains for this rather complicated interplay 
of forming and stabilizing extended textures and aligning domains. Apparently, the use 
of this growth law for estimating domain sizes seems justified throughout the whole time 
period relevant for DCC signatures. 

It is interesting to compare these results with the plain linear cr-model. For that purpose 
we omit in (0) the current-current coupling and follow the evolution without putting a 
constraint on the winding number B. Fig.|] shows two events evolving from identical initial 
configurations on a 100x100 grid with vacuum boundary conditions. After about 200 relax- 
ation time units the unconstrained linear sigma model has lost all particles and antiparticles 
{D — > 0), and the configuration tends towards global alignment in 3-direction. On the 
other hand, the full model with fourth-order Skyrme term included and winding number 
constrained to -B = 0, finally stabilizes solitons with D of the order of 100 particles plus 
antiparticles. Comparing the growth rates of the angular correlation lengths Rjj for both 
evolutions it may be seen that there is no dramatic difference in the overall growth exponent. 
Still, it is evident that especially during the roll-down phase the emerging extended stabiliz- 
ing textures do present an obstacle against the free and unconstrained ordering process such 
that the average size of the ordered domains in the full model is noticeably smaller than in 
the plain linear a-model evolution. 
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FIG. 3. Comparison of particle numbers D, angular correlation lengths Ro^ and the average 
bag fields (<1>) as functions of time t after a sudden quench, for evolutions under the full model with 
constraint i? = (full lines), and the pure linear a- model (no fourth-order term, no constraint on 
B (dashed lines); (100x100 grid, l = 2,\ = 2). 



At this point we may add a few remarks concerning the relaxation constant r which in 
the TDGL-equation (|18D we have simply used as a time unit, because no other propagating 
terms are present. Of course, we could use (^) instead of (|^) for simulating evolutions of the 
classical field. Keeping the propagation velocity at c = 1 the time unit then is set by the 
lattice constant. A trivial consequence is that all features which reflect the relaxation, like 
growth rates or saturation times, are scaled with powers of r. E.g., for sufficiently small r 
the roll-down time Ird scales like 1/r. It turns out that for r ~ 0.1 the evolutions resulting 
from (^ and (|^) essentially coincide. For r ~ 1 the propagating waves start showing up 
in the field configurations and cause oscillatory features in (lattice-) averaged quantities 
without strongly affecting the late-time evolution. For r > 2 oscillations dominate and due 
to the nonlinear terms stabilization becomes difficult. So, to be on the safe side for use of 
the TDGL equation (^, r should be less than 1. Looking once more at fig.^, tRr, ~ 10 (in 
units of r) implies that for r ~ 0.1 the roll-down is completed after travelling waves have 
passed about 100 lattice units. 

V. COOLING TIMES 

The assumption of a sudden quench which underlies the ordering features shown in the 
previous chapter represents that limit for the actual cooling process where the cooling time 
Tc (during which the temperature drops from T /T^ ~ 1 to T = 0) is much smaller than 
the relaxation time which governs the evolution of the field configurations. This instanta- 
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neous cooling, where the hot initial configuration for t > is immediately exposed to the 
zero-temperature T = effective potential defines the most off-equilibrium scenario which 
we can study within this model. To weaken this rather drastic assumption we replace it 
with an appropriate time dependence of the temperature T{t). Then also the temperature 
dependence p{T) in (pj| ) has to be specified. 

The function P{T) in principle is determined from loop contributions and their renor- 
malization group resummation in finite-temperature field theory. Perturbative evaluations 
are restricted to the vicinity of T = 0; the complete form of P{T) depends on the specific 
model and is not accurately known. Its essential features, however, are a slow decrease from 
its T = value /^ = 1 for T > 0, followed by a steeper decrease /^ — ^ as T approaches 
Tc- For the simulations these essentials are sufficient; the specific form (generally in form of 
logarithms) is not important. A very convenient ansatz which comprises these basic features 
in the region < T < T^ is 

P{T) = (1 - {T/TX)f (24) 

with a ~ 4, and p{T) = for T > T^. 

A slow cooling scenario can be realized by externally imposing the dependence T{t) of 
temperature on time. For a linear quench we take 

T(t)/T, = l-t/r, (25) 

with quench time Tc-, and T = for t > r^. A crucial feature of the ordering process is that 
the early phase prior to the onset of the roll-down is dominated by angular aligning driven 
by the gradient terms in the a-model part of (0) while the effective potential for $ is of only 
minor importance. As a consequence, the angular ordering up to the onset of the roll-down 
proceeds quite independently of the actual cooling rate, with ($) <^ 1 in any case. This is 
fortunate because it renders knowledge of the precise form of T(t) and P{T) unimportant 
and we can safely adopt convenient ansaetze like ( P^ and (P5|). 

For a very slow quench Tc 3> r, after this first angular aligning phase, the roll-down then 
only proceeds up to the momentaneous value of /^, (<l')^ < p(T(t)). Specifically, the stable 
textures are being created with profiles which correspond to an instantaneous scale 

m = if-\T{t)) (26) 

as can be seen from rescaling the modulus of the order parameter to $ = <l'/^^(T(t)) in (|I1|). 
With p{T{t)) < 1 the nascent stable objects in the slow quench are much larger than for 
the sudden quench, i.e. fewer of them survive as separate isolated textures. Subsequently, 
further ordering proceeds essentially as equilibrium evolution where the stable objects shrink 
according to £(t) in (^) until they finally assume their T = size with accompanying 
rearrangement of surrounding aligned domains. 

As intermediate case we consider evolutions where the cooling time is comparable with 
the roll-down time Ird in (0); conveniently again as linear quench with Tc ~ 10 for £ = 2 
and A = 2. In the very early phase the evolutions closely follow the slow quench with 
very small ($) and steeper decrease of the particle number D as compared to the sudden 
quench, but as t approaches Tc the roll-down proceeds rapidly towards the sudden quench 
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FIG. 4. Energies E, particle numbers D, angular correlation lengths Rd, and average bag fields 
($) for three events (lOOx 100 grid, A = 2, ^ = 2) with linear quench (24) and three different cooling 
times: Tc = (sudden quench), Tc = 10 (intermediate), Tc = 100 (very slow quench). The dashed 
line indicates the growth law Rd ~ t^'^. 

value ($) ~ 1, with corresponding pronounced loss in energy E and typical increase of Rd- 
For t > tRD the number of textures with T = profiles stabilizes at values intermediate 
between sudden and very slow quench and the angular ordering of aligned vacuum domains 
proceeds along the classical t°'^-law. 

Fig.| shows three typical evolutions with i = 2 and A = 2 and total winding number 
constrained to -B = 0, starting from the same initial configuration, for r^ = 100 (very 
slow), Tc = 10 (intermediate), and Tc = (sudden) quench. Fig.|^ presents bag contours at 
$ = 0.4, 0.6, 0.8 of the stable textures with varying individual winding numbers (summing 
up to -B = 0) which have been established after t ~ 100 for these three different cooling 
rates. It is quite evident from fig.|^ that the rate of increase of the correlation lengths Rd 
before and after the roll-down (t ~ 10 — 20) is unaffected by the cooling rate. Differences in 
the angular ordering, however, show up in the evolving particle numbers which reflect the 
changing effective potential through the size of nascent textures. 

For these slow quenches which proceed with finite cooling rate evolutions are subject 
to fiuctuational noise given by the last term in (|18]). Its variance depends on the physical 
nature of the dissipative $ term in (18), i.e., to which extent it originates purely from 



particle emission, or from kinematical expansion of the spatial volume, or from a heat bath 
of eliminated high-frequency modes. In any case, its strength according to the fiuctuation- 
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FIG. 5. Bag contours (drawn at <& = 0.8,0.6,0.4 from outside to inside) of isolated baryons 
remaining at time t = 100, for the 3 events shown in fig.^ which differ by their coohng times 
Tc = 0, 10, 100. By the time t = 100 the temperature has arrived at T = in all 3 cases, and 
the baryons approach their T = sizes and quasistable shapes. The smallest ones carry one 
unit of winding number. The total winding number is zero in all three cases. The remaining 
particle numbers (as can be seen from fig.§) are D = 95, 65, 24 for the cooling times t^ = 0, 10, 100, 
respectively. 



dissipation theorem decreases linearly with temperature T. Its general influence on the 
evolving configurations is very transparent: for finite T it will tend to disturb the ordering 
process at each time step. This results in a slight decrease of the slopes of those observables 
which reflect the ordering, specifically the correlation length Rd, the average ($), the particle 
numbers D, the total energy E. On the other hand, fluctuations support (Brownian) motion 
of nascent textures which combine more easily into larger structures with increased individual 
winding numbers. So, evolutions subject to noise finally produce increased values of D 
distributed into larger clusters of positive or negative winding number. Naturally, only very 
slow quenches are noticeably affected by the noise. For definiteness, the evolutions shown in 
figs.^ are obtained from (|T^ with the last term $, omitted completely. Effectively, inclusion 
of noise produces evolutions which resemble those for quenches with slightly reduced quench 
time. 

As we have seen, the physically most relevant characteristic of the ordering evolution is 
the roll-down time where the bag field $ approaches within aligned (but still mostly disori- 
ented) domains the momentaneous value of f(T) and isolated textures are being established 
near their borders. This important feature is most clearly brought forth for r^ ^ therefore 
in the following we preferably will consider ordering processes after instantaneous cooling. 



VI. EXPLICIT SYMMETRY BREAKING 

Explicit symmetry breaking of the global 0(3)-symmetry is of crucial importance for the 
ordering process. Its weakest form is imposed through boundary conditions where a specific 
vacuum orientation, e.g. $ = (0,0, 1), along the lattice boundary is enforced. For physical 
processes this, apparently, is also a most natural condition. With increasing time after a 
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sudden quench the influence of this surrounding vacuum will grow deeper and deeper into 
the randomly ordering interior and eventually reorient aligned domains along the chosen 
vacuum direction. Clearly, the strength of this aligning effect depends on the finite size L of 
the lattice. For L of the order of 100 it does not change the essential features of the evolution. 
There is only one ordering parameter, the absolute value of the averaged 'direction' 



L-l 



W= Y.H^^j)/L' 



(27) 



i,j=0 



which is sensitively affected. For arbitrary periodic boundary conditions this observable 
remains very small with erratic variations for times long after the roll-down, while for the 
symmetry-breaking vacuum boundary conditions, it continually grows in accordance with 
the growth of the angular correlation length. It is saturated by the averaged 3-component 
($3) which carries all of the total (<&) strength. So, as expected, the oriented domains 
growing in the vicinity of the boundary contribute coherently to the magnetization and 
align in vacuum direction. 




0.1 1 10 

FIG. 6. Angular correlation lengths Rd^ average bag fields {^), and the averaged 3-coniponent 
of the field unit vector (^3), for three events (100x100 grid, A = 2,£ = 2) with different strengths 
of explicit symmetry breaking, m^ = (full lines), m^ = 0.01 (dashed lines), m^ = 0.1 (dotted 
lines). (Note that the average bag fields (<I>) for m^ = and m\ = 0.01 almost coincide.) 

Explicit symmetry breaking through a nonvanishing m^-term in ([T^) exerts much 
stronger influence on the ordering process than imposing symmetry-breaking boundary con- 
ditions only. We compare in fig.§ evolutions under vacuum boundary conditions for values 
of m^ = 0.01, 0.1, and ml = 0. This appears as a reasonable range for pions with mass 
of about 10% of the a-mass. For m'^ = 0.01 the roll-down of ($) is almost unaffected as 



compared to m^ = 0, while m^ = 0.1 accelerates the roll-down by a factor of two (with a 
correspondingly larger number of stable textures created). 

Both evolutions with m"^ ^ show the exponential approach of ($3) to saturation 
expected for the off-critical quench, in contrast to the unbroken case. (For infinite systems 
the exponential approach to saturation is governed by exp{—2m'^t) |1^). The saturation 
value is smaller than one due to the presence of stable textures. In both cases {^) is 
completely saturated by ($3). It may also be noted from fig.|^ that for m^ = 0.1 saturation 
of ($3) is achieved long before the roll-down of ($) is completed, while for m'^ = 0.01 the 
order is reversed. 

It should be noted that the angular correlation function C{R) in the off-critical system 
approaches ($)^ for R ^ 1. This shows up as an abrupt increase of the half-maximum 
distance Rd, as soon as the average ($) has grown beyond -\/0.5 ~ 70% of its saturation 
value. Before that sudden increase, Rd already starts to deviate from the growth law of the 
critical system, because it reflects also the increasing global alignment in 3-direction. So, in 
the off-critical system, it appears more appropriate to consider only the angular correlations 
of the transverse {i = 1,2) components as a measure for the size of disoriented domains. 

Additional breaking of the 0(3)-symmetry occurs through an initial bias in the starting 
configurations. It originates from the asymmetry introduced into the initial T > Tc ensemble 
through the m^$3 symmetry-breaking term in (0). It will not alter the exponent of the 
exponential approach to saturation, but it will reduce the factor in front of the exponential. 
This causes another considerable reduction of saturation times. 

VII. PION EMISSION AS DCC SIGNATURE 

As we discussed in sect. 2 signatures of the evolving disoriented domains may be carried 
in the fractions /(vtj) of pion multiplicities of kind i emitted at late time during the ordering 
process. We will investigate these signatures for an ensemble of events evolving after a 
sudden quench. We compare results for the plain linear a-model, (i.e. omitting the four- 
derivative current-current coupling in (pjlf ), without constraint on winding number B), with 



the corresponding evolutions for the full lagrangian (p!^, with B constrained to i? = 0, 
which lead to stabilized solitons floating in the surrounding vacuum. For the parameters in 
the energy functional we choose i = 2 and A = 2, such that the emerging solitons develop 
deep bags and extend over several lattice cells. With this choice the mass of the (T=0) 
(j-fluctuations according to (|T^) is rria- = 1. 



Let us at first look at the rates e^ ([TlD of energy carried away by the field components 
i = 1,2,3 for the pure linear cr-model with no explicit symmetry breaking, ?7i,r = 0. The 
direction of spontaneous symmetry breaking then is fixed by the surrounding vacuum alone, 
which is imposed through the condition that the chiral field is kept at $ = (0, 0, 1) at all 
times on the boundary of the lattice. For these events the roll-down phase where the average 
bag field ($) grows from a small value to saturation ($) ^ 1 covers the time intervall from 
t ~ 1 to t ~ 10 in units of the relaxation time r. During this time interval the correlation 
lengths Rd grow from values near Rd ~ 1 to Rd ~ 10. These growth rates are quite 
independent of the lattice size L, such that the average number of disordered domains by 
the end of the roll-down phase is of the order of (L^/(7ri?|,). 
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FIG. 7. Angular correlation lengths Rd, average bag fields (<1>), and the rates e^ of energy 
loss ( [Tl|) carried by each of the three components i = 1,2,3 of the chiral field, as functions of 
time t after a sudden quench, on a 30x30 (left) and an 80x80 (right) grid (plain linear a- model, 
no constraint on B). 



In figs.|^ we show for some arbitrarily selected individual event the rates e^ of emitted 
energy for two lattice sizes, L = 30 and L = 80. Due to the symmetry breaking effect of 
the surrounding vacuum, emission into the i = 3 a-direction dominates the pionic i = 1,2 
emission rates during most of the roll-down phase by one to two orders of magnitude. 
Naturally, this feature is more pronounced in the smaller lattice. Only after completion of 
the roll down, when the absolute emission rates are getting very small, the average rates 
become again comparable for i = 1,2,3. While €3 is essentially monotonously decreasing 
during the whole ordering process, the pionic emission occurs in two well-separated pulses. 
An early, short pulse with high intensity is emitted for times t < 1 and it is characterized 
by equipartition ei = €2- This is as expected because at those times the angular correlation 
lengths Rd are near or less than 1. The second pulse appears during the second half of the 
roll-down phase; its maximum is at least one order of magnitude smaller than that of the 
first pulse and it is reached near ($) ^ 0.8. For the larger lattice also this second pulse is 
characterized by equipartition for both isospin components. However, for the smaller lattice, 
a significant difference between ei and €2 can be observed. It develops with the onset of the 
roll down and extends over the whole second pulse of pion emission and thus constitutes a 
possible DCC signature. Unfortunately, it is just in the small lattice where the pionic signal 
is strongly dominated by a-emission. (Counting particles instead of energy, this is, however, 
compensated by the mass ratio). 

It is instructive to analyse the two pion pulses in terms of the TDGL-equation (|18|) 
which relates (in the pure linear a-model with m^ = 0) $i to A$j — (A/£^) ($^ — /^) $j. 
As expected, the first pulse is mainly due to the angular ordering effectuated by A<l>j, 
while the second pulse is dominated by the potential term. So, the simple approximation 
e, ~ (($^ — p) $f) provides the essential features of the second maximum, (the 'DCC 
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FIG. 8. Probability distributions P{f) for the abundance ratios fi{t) (|13D of DCC pions ex- 
tracted from an ensemble of 3000 events at roll-down time i_R_D. For comparison, suitable mixtures 
of Pu{f) selected from the approximation ( |A5| ) are shown by the dashed lines. The events evolve 
after a sudden quench on a 30x30 grid, with no explicit symmetry breaking, according to the plain 
linear cr-model without constraint (left), and for the full model with B = Q (right). 
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FIG. 9. The same as in figj^, but for a 60x60 grid. 

pions'). We may even go further and replace $^ by its lattice average ($^), such that the 
abundance ratios (0) for DCC pions are approximated by 



f{^^ 



(*^ 



m) + ($1) 



(28) 



In this approximation we recover the commonly used prescription [Q to take the strength of 
pion emission proportional to the square of the classical field components. (Evidently, this 
works only for ratios, and only for pions in the second peak). Based on this approximation, 
we may derive simple expressions for the probabilities Pu{f), to observe in a large ensemble 
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of configurations with v disoriented domains the ratios /(vtj) of DCC pion multiphcities (see 
Appendix). 





FIG. 10. Average bag fields (<l>), averaged 3-coniponents of tlie field unit vector (^3), and the 
rates Ej of energy loss (11) carried by each of the three components i = 1,2,3 of the chiral field, 
as functions of time t after a very slow quench {jc = 100), on a 30x30 (left) and an 80x80 (right) 
grid (plain linear cj-model, no constraint on B). 

As discussed in the previous chapter, addition of exphcit symmetry breaking through a 
small pion mass like m^ = 0.1 speeds up the roll-down phase and the aligning of $ into the 
3-direction and thus shifts the DCC pions peak to slightly smaller times without affecting 
the ratios. However, the important difference is that for 771^^ 7^ the emission rate for the 
3-component is strongly enhanced relative to the DCC pions, such that also for a large 
lattice {L = 80) during and beyond the whole roll-down phase es dominates the DCC pions 
by more than an order of magnitude (for m.,^ = 0.1). (For the small lattice {L = 30) €3 
dominates by two orders of magnitude). 

The fourth-order Skyrme term plays a minor role during these early stages of the ordering 
process. Its size stabilizing effect for the emerging solitons affects the late stages after 
completion of the roll down. Still, as we have discussed at the end of sect. 4, already the 
emergence of stable extended structures during the roll-down phase presents an obstacle 
against the growth of large aligned domains and thus it counteracts large differences in 
the ratios /i and /2 of DCC pions. This becomes apparent by comparing the probability 
distributions P{f) obtained from a large ensemble of 3000 events on a 30x30 lattice (fig.^, 
for the plain unconstrained linear a-model, and the full model (|1^ with constraint B = 0. 
We sample into 30 bins the ratios /(vti) extracted from 3000 evolutions at the time when ($) 
grows beyond 0.8, which approximately marks the peak in the DCC-pions pulse. Comparing 
with linearly mixed Py{f) taken from ( [A5| ) for two neighbouring values of v, we may conclude 
that the presence of stable textures increases the average number of disoriented domains 
present at that point in time (for the fixed lattice size L = 30) from 2-3 to 4-5. As we have 
discussed above, this average number of u increases with L^. Fig.|| shows the corresponding 
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distributions for 3000 events on a L = 60 lattice. The distribution obtained from the hnear 
cr-model is very nicely reproduced by Pu{f) with u = 13, while the full and constrained to 
B = evolutions are compatible with i> = 27, which again shows the strong obstruction 
of domain growth through the creation of stable textures. Still, it should be remembered 
that the Pu{f) of ( |A^ ) which we use for comparison are derived under the rather poor 
approximation (^) and the assumption that all domains are of equal size. 

To illustrate the influence of the cooling rate for slow quenches figs.|T^ show two evolutions 
with cooling time Tc = 100 as discussed in sect. 5. With this choice the cooling time is large as 
compared to the roll-down time after a sudden quench. The decisive feature of the evolution 
on a small 30x30 lattice is that the saturation of ($3) is almost completed before the begin 
of the roll-down of ($). Consequently, the DCC pion pulse which still occurs during the late 
roll-down phase, is suppressed by almost four orders of magnitude as compared to cr-emission 
63. On this small lattice the disoriented domains show up in the pion emission ratios, while 
on a larger (80x80) grid, the saturation of ($3) is delayed, the height of the DCC pion pulse 
increased, but it is characterized by equipartition between ei and €2- So, also in this respect 
the sudden quench provides the better conditions for a DCC signal. 

VIII. CONCLUSION 

In this paper we have investigated the phase ordering of a 2D-0(3) model under dissi- 
pative dynamics. The peculiarity of the underlying lagrangian is that it stabilizes extended 
textures through a fourth-order current-current coupling. This is introduced to get an idea 
how emerging topological solitons with stable sizes interact with the growth of aligned do- 
mains. The winding number can be kept fixed to simulate in this model the analogy of baryon 
number conservation which is an important feature of models of hadronic systems. Two pa- 
rameters characterize the model: One fixes the spatial size scale of the stabilized 'baryons', 
the other determines the depth of the bag profile, into which the winding (baryon-) density 
is confined. Shallow bags correspond basically to the nonlinear version of the model, deep 
bags realize the unbroken chiral symmetry in the interior of baryons, as it would be expected 
in hadronic systems. 

Apart from studying the influence of different ways of symmetry breaking, and the effects 
of different quench rates on the ordering observables, we are especially interested in the 
possibility to find DCC signatures in the emitted particle radiation. All results are based on 
the assumption that the relaxation proceeds on a sufficiently slow time scale, such that it 
provides an adiabatically changing source for the emitted radiation. Therefore its motion is 
treated classically, and the amount of radiated energy is simply obtained from the classical 
energy loss of the source. Quantum field-theoretic amplitudes enter only into the relaxation 
constant r which, however, has not been determined, but simply used as time unit for the 
ordering dynamics. 

From the results obtained in this framework we may draw the following conclusions: The 
classical AUen-Cahn growth law Rd ~ t^'^ for disoriented domains is not seriously disturbed 
by the stabilizing terms in the lagrangian and the emerging stable textures, although there 
are characteristic deviations from it during the roll-down phase where the separation into 
domains of false vacua and deep bags mainly takes place. The factor in front of the t^'^ power 
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law, however, is decreased as compared to the pure unconstrained hnear a-model, such that 
the average size of disoriented domains after completion of the roll-down is reduced. The 
roll-down time varies with the two parameters of the model essentially as expected from 
simple scaling arguments. It represents the crucial characteristic of the ordering process, it 
is not very sensitive to small explicit symmetry breaking and varies with the quench rate in 
the expected way. Explicit symmetry breaking strongly affects the saturation time of the 
average orientation, the factor in front of the exponential approach to saturation depends 
severely on the (random) initial conditions of each individual event. The number of baryons 
(plus antibaryons) finally surviving the ordering process is determined by the time of the 
onset of the roll down. The angular correlation length at that time measures the average 
distance of regions with high winding density which develop into separate stable structures 
with integer winding number during and after the roll down phase. Therefore slow quenches 
lead to low baryon density. 

The intensities of emitted radiation in our approach is determined by the square of the 
time-derivative of the classical field, averaged over the spatial volume (or area), and not by 
the square of the classical field itself. This allows to identify two well-separated radiation 
pulses for the transverse ('pionic') components of the chiral field. The later one, emitted 
towards the end of the roll-down phase carries the signature of the disoriented domains 
established by that time, so it constitutes the DCC signal. With the growth law of the 
angular correlations fixed, the average size of disoriented domains at the time of completion 
of the roll down is also fixed. Therefore the average number of disoriented domains on an 
L X L lattice which contribute to the DCC signal increases like L^, and deviations from 
equipartition in the intensities carried by the different isospin directions quickly vanish with 
increasing L. On the other hand, for small L, the pion signal is buried under a large pulse 
of emission in longitudinal ((T-)direction, because the explicit symmetry breaking effect of 
the surrounding true vacuum is strong for small lattice sizes. 

The broadening of the probability distributions for the pionic intensity ratios consti- 
tutes the DCC signature. We have extracted it from large ensembles of events under the 
assumption that these pions, emitted at a specific point in time (or averaged over a short 
time interval near the roll-down time) can be observed separately from the bulk of other 
pions and cr's emitted. We find that in those cases where a noticeable broadening occurs, 
the relative amount of energy carried away in the DCC pulse, is suppressed by orders of 
magnitude. Of course, quantitative statements require knowledge of the parameters in the 
model, cooling rates, physical size of the lattice as compared to correlation length, baryon 
size. In the present model, all these conditions remain more or less arbitrary, and we have 
chosen them at will to illustrate the main features. However, for application to the ordering 
dynamics of a hadronic chiral 0(4)-field, we have quite specific ideas about the relevant pa- 
rameters in the effective actions. The exponent in the classical Allen-Cahn scaling law is the 
same for (2+l)D-0(3) and (3+l)D-0(4) models; the homotopy groups 7i2{S^) and 7r3(5'^) 
are isomorph, both characterized by integer winding numbers. So, apart from differences 
in dimension-dependent cooling rates and temperature-dependence of the condensate /(T), 
we do not expect dramatic qualitative differences in the ordering features, and we hope that 
extending the present considerations to the (3-|-l)D-0(4) model will serve to sharpen our 
expectations with respect to the observability of DCC signals in hadronic processes. 
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APPENDIX A: PROBABILITY DISTRIBUTIONS 

Consider a configuration where the classical field is uniformly aligned into some direc- 
tion, thus forming one single disoriented domain, with the field components {i = 1, 2, 3) 
parametrized as 

^j=i,2,3 = $(sin 9 cos (j), sin 9 sin 0, cos 9) . (Al) 

Then, within the approximation (pSf), we find for the fraction of pions emitted with isospin 
component i = 1, relative to the number of all pions {i = 1, 2) 

/(7ri) = cos2 0. (A2) 

In an ensemble of events where all orientations of <l>i_2 are equally probable the ensemble 
average (/) of /(vti) is, of course, (/) = 1/2. The probability P{f) to find in one event of 
that ensemble the fraction / (for pions of either kind i) then is obtained from 

U) = ^ r fdct^= f fPU)df (A3) 

zvr JO Jo 



as 



P{f) = I J ■ (A4) 

(Note the difference to the standard chiral 0(4)-model, where P{f) = l/{2y/f).) For u such 
domains of size 1/z/ the probability Pu{f) to find the fraction / in a given event then is 

P,(/) = £ ... £ 5^f- \{h + - + fu)) P{fi)...P{f.) dh...dU (A5) 



with P{f) given by ( |A^) . Some of these functions, for values from z/ = 1 to i^ = 10, are 
shown in fig.|ll|. For increasing v they are, of course, more and more sharply peaked around 

/ = 1/2. 
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FIG. 11. The probability distributions Pu{f) for observing the fraction / of pions with a given 
isospin component i = 1 or z = 2 in an ensemble of configurations with u disoriented domains of 



size 1/z/, according to (A5). 
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